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We have developed a theory of polymer entanglement using an extended Cahn-Hilliard func- 
tional, with two extra terms. One is a nonlocal attractive term, operating over mesoscales, which 
is interpreted as giving rise to entanglement, and the other a local repulsive term indicative of ex- 
cluded volume interactions. This functional can be derived using notions from gauge theory. We 
go beyond the Gaussian approximation, to the one-loop level, to show that the system exhibits a 
crossover to a state of entanglement as the average chain length between points of entanglement 
decreases. This crossover is marked by critical slowing down, as the effective diffusion constant goes 
to zero. We have also computed the tensile modulus of the system, and we find a corresponding 
crossover to a regime of high modulus. The single parameter in our theory is obtained by fitting 
to available experimental data on polystyrene melts of various chain lengths. Extrapolation of this 
fit yields a model for the cross-over to entanglement. The need for additional experiments detailing 
the cross-over to the entangled state is pointed out. 

PACS: 61.41.+e, 83.10.Nn 



I. INTRODUCTION 

While it has long been known that entanglement in 
homopolymers has an important effect on its strength, a 
thoroughly satisfactory theory of polymer entanglement 
is still a topic of current research, ffhe classic exper- 
imental work of Moore and WatsonEI showed that the 
bulk modulus of cross-linked natural rubbers depends in- 
versely on the average chain length (N c ) between cross- 
links in the system, and that end corrections become 
negligible as the total average molecular weight [i gets 
very large. They pointed out the analogy between chem- 
ical cross-linking and physical entanglement. Thus their 
work applies in a qualitative sense to entangled systems 
as well. Their, work extended the earlier pioneering work 
of Flory et afl'Q. 

Edwards developed the tube theory of the effect of en- 
tanglement on elastic moduli of homopolymers using de- 
Genne's idea of reptatiora. This theory showed how en- 
tanglement enhances the tensile modulus of a homopoly- 
mer. He also developed a more detailed model of en.- 
tangled ring polymers using notions from knot theoryQ. 
The basic idea behind this theory is an analogy between 
certain mathematical invariants describing intertwined 
loops and magnetic fields induced in wires by current- 
carrying loops. Prager and FrischB workcd-on this notion 
as well, as did Koniaris and MuthukumarO. 

More recently, interest has turned towards computer 
simulations of polymer networks, involving various levels 
of molecular detail, to understand the effect of entangle- 
ment on the strength of homopolymers. As exajnples, we 
mention the work of TermoniaQ and Bicerancu, who use 
phenomenological models of polymer networks to study 
their viscoelastic properties. Comparison with experi- 
mental data shows a varying degree of success, depending 
on the particular system studied. Holtzl et alcl use the 



more basic fluctuating bond theory to model a network 
of polyethylene strands to show that entanglement leads 
to non-afnne displacements under large tensile strains. 

In an earlier paper ,E3 we developed a gauge theory of 
self-assembly, and utilized renormalization group ideas to 
study the onset of self-assembly in diblock copolymers. 
In this paper, we shall pursue a similar continuum ap- 
proach to understand entanglement. 

Intuitively, one can see that entanglement could be de- 
scribed by .assuming two extra terms in the Cahn-Hilliard 
functionaljli one of which is a nonlocal attractive term, 
which gives rise to entanglement and the other, a soft- 
core local repulsive term which arises from the fact that 
the strands cannot cut across each other. We connect the 
parameters which appear in our theory to the underlying 
chain parameters with a simple model. We have shown 
(see Appendix A) how such a functional can be derived 
naturally using notions from gauge theory. 

The results derived from our continuum formulation 
will be seen to be reminiscent of the. chain-theory ap- 
proaches of Kavassalis and NoolandO for flexible poly- 
mer networks, and that of Kroy and FreyLJ for semi- 
flexible networks. They utilized a mean-field approach to 
locate the transition to the state of entanglement. Our 
theory is alsa.somewhat similar to the paper of Castillo 
and GeldartEfl, who use a 4? field theory, coupled to the 
replica trick £in the mean field approximation) of Deam 
and EdwardsEj to study the vulcanization transition. 

We shall utilize a field theoretic approach and go be- 
yond the Gaussian approximation in this paper to show 
that the onset of the state of entanglement is a crossover 
phenomenon, rather than a pure phase transition, in that 
the effective diffusion constant goes to zero at the tran- 
sition point, but the correlation length and the structure 
factor do not diverge. A physical reason that fluctuations 
become important near the onset of the state of entan- 
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glement is that the average chain length between points 
of entanglement gets smaller, while in a vulcanized poly- 
mer, the cross-links make the system increasingly stiff. 
This underscores a difference between vulcanization and 
entanglement- An entangled network of polymers is more 
dynamic than a vulcanized network. The mean field ap- 
proximation is expected to be correctu for the vulcan- 
ization transition. 

We have also computed the tensile modulus of the sys- 
tem. Corresponding to the critical slowing down dis- 
cussed above, we find a crossover in the modulus to a 
regime of high values. Fits to available experimental 
data shows how the single parameter in our theory can 
be parameterized in terms of the molecular weight of the 
system. Extrapolation of the parameterization we have 
provided in this paper displays the cross-over to the en- 
tangled state. The need for further experiments detailing 
this cross-over is pointed out. 



II. A FIELD THEORY OF ENTANGLEMENT 

The continuum, mesoscale approach adopted in this 
paper assumes that we have performed some spatial av- 
eraging of our polymeric system, so that the order pa- 
rameter is the local concentration of the polymers. Our 
mesoscopic theory of entanglement in polymers is based 
on the intuitive notion that physical entanglement can be 
captured by a non-local attraction between the polymers 
which causes them to remain in proximity. There must be 
a balancing, repulsive local energy term which says that 
the polymers cannot cut across each other. The start- 
ing point of our mesoscale theory is an internal energy 
functional which is quadratic in the gradient of the local 
number concentration. For the moment, we will consider 
isolated systems, so that the quantity that is conserved 
is the internal energyiiH. We will shortly consider entropy 
effects as well. Consider the following form for the energy 
functional: 



/3U =P « (c(s))d J 
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Pu (c(s)) = (I) 



dc(s) dc(s) 
dsi dsi 



(1) 



(2) 



(3) 



where repeated indices are summed over, s is a dimen- 
sionless co-ordinate variable, k is Boltzmann's constant, 
T is the temperature, and c is the number concentration 
of the specie. The local concentration c is normalized 
by dividing by some characteristic inverse volume. The 
constant g is analogous to a dimensionlcss diffusion con- 
stant. Such energy functionals have been considered over 
many years as contributing to the total internal energy of 



both unary and binary mixtures0 We will use this form 
as our starting point to suggest a more complete energy 
functional: 



(3U cS = (3U + 



J d 3 s c(s)c(s) 



c(s') (4) 



where a 2 , 7, 6 are positive constants. The local repul- 
sive term is indicative of the fact that polymers cannot 
cut across each other. This is in effect a soft-core repul- 
sion term, and the softness arises because we are study- 
ing a homopolymer network at a mesoscale, where poly- 
mers may pass by each other, without actually cutting 
across each other. The nonlocal attractive term gives 
rise to entanglement, as it causes portions of the net- 
work within the screening distance 1/5 to be attracted 
to each other. Equation ^ is the basic statement of our 
theory. Note that the two terms we just discussed have 
signs opposite those of corresponding terms in theories 
of self-assemblyll3. The non-local attractive term invokes 
the notion that entanglement of polymers must lead to 
knotty configurationsOj. This attractive term can be seen 
intuitively to lead to the notion of trapping, and as such 
is slightly similar to treatments of entanglement in the 
literature which use Arrhenius-like rate theories to pro- 
vide for escape of polymers from local entangled arrange- 
ments. 

In what follows, we shall set 7 = a 4 , and S 2 — \j2a, 
with a 2 = g 2 /2. A strong motivation for this choice of 
parameters is provided in Appendix A, where we use no- 
tions from gauge theory to derive Eqn^, with the param- 
eters having the forms given above. Another explanation 
for such a choice is as follows. With our choices for the 
parameters, U e fi in momentum space may be written as: 



d 3 k 
(2^)3 



*(fc)[-\/2afc 2 - V2ak 2 /(1 + V2k 2 /a)]c(k) 



(5) 



where the carats indicate a Fourier transform. Thus we 
see that the choices made for the parameters are equiv- 
alent to generalizing the diffusion constant g = y2a — > 
\/2a[l + 1/(1 + y/2k 2 /a)], i.e., a non-local diffusion con- 
stant is obtained. If we now extremize the functional, the 
Euler-Lagrange equations can be written in conservative 
form as: 



V • I(s)= 
I(s)= V 



d 3 k 



■ exp 



(ik- s)[l + + V2k 2 /a)]c{k) (6) 



where I(s) can be interpreted in the conventional man- 
ner as a mass current. The divergence-free nature of 
this current makes it clear that with our choice of pa- 
rameters, our internal energy functional preserves num- 
ber conservation. This is quite appropriate, since the 
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internal energy is the quantity which is conserved for iso- 
lated systems. For an arbitrary choice of parameters, 
the Euler-Lagrange equations have the form: V • I'(s) = 
Source/Sink Terms, indicating that number conservation 
can be a problem. 

While our choice of parameters may appear to be 
overly restrictive, it turns out to be sufficiently rich to 
provide a description of the onset of entanglement in 
polymers. We will not explore more general sets of pa- 
rameters in this paper. 

Before we can compare our theory with experimental 
data, we need to consider the fact that our system is not 
really isolated, and may be in contact with an energy 
reservoir, perhaps as it is being acted on by mechanical 
forces in a stress experiment. For a system in contact 
with an energy reservoir, the-quantity that is conserved 
is the Hclmholtz free energyO A — U c s — ST, where S 
is the entropy of the system. The entropy of our .system 
will be approximated in the conventional mannertj: 

- | = J d 3 S c(s) ln[c(s)] (7) 

Note that we have ignored a term in the above ex- 
pression which is linear in the normalized concentration 
field. This term can be absorbed into the definition of 
the usual Lagrange multiplier constraint for number con- 
servation. This additional constraint is necessary, over 
and above the considerations which led to Eqn.||, because 
we are now considering the Helmholtz free energy rather 
than just the internal energy. In the homogeneous mean 
field approximation, the chemical potential can be easily 
shown to be zero. We shall utilize this approximation to 
facilitate computations. This entropy term provides the 
free energy a single minimum. To further ease computa- 
tions, we shall expand cln(c) in a power series about the 
characteristic inverse volume ^ _3 (= 1) in our dimension- 
less units), retaining terms up to fourth order: 

(l + C )ln(l + C )« C +^-^ + ^ (8) 

This expansion yields one minimum, just as the exact 
expression for the entropy. Consequently, we do not ex- 
pect this system to display a phase transition, but rather 
a crossover from an un-ent angled state to a state of en- 
tanglement. Finally, we note that in our present theory, 
entropy yields the crucial nonlinear terms which will de- 
scribe the crossover to a state of entanglement, in con- 
trast to our gauge theory of self-assemblytB where en- 
tropy did not play a dominant role. 

We define the two-point Green's function as usual via 
S{x,x') = Limj^ (5 2 /SJ(x)SJ{x'))Q[J], where Q[J] = 
J VcB{l + c)exp[-j3(U eS - ST) - J d 3 sJ(s)c(s)}, where 
6(1 + c) is a step function that indicates a restriction 
to physically acceptable values of the concentration. In 
practice, we shall be restricting our attention to small 
deviations of c from its average, so that the step func- 
tion is implicitly accounted for during calculations. In 
the quadratic approximation, the structure factor is: 



0[ ' l + a'k 2 + a' 2 k 2 /(l + 2k 2 /a>) ( ' 

where a' — y/2a. Equation |^ displays a peak at the 
origin, as one might expect from the fact that entan- 
glement creates blobs which are distributed at random 
within the system. The width of the peak indicates an 
inverse of the correlation length between the blobs. With 
this physical interpretation, yfo. is a measure of the dis- 
tance between concentration fluctuations (i.e. between 
points of entanglement). The decay of So(k) is affected 
by the value of a. As a decreases, the structure factor 
looks more diffuse. Thus, a decrease in a signifies a shift 
to a state of higher entanglement, as the concentration 
of entanglement points increases. 

Our results can be understood compactly in terms of 
the parameter a, or equivalently, g, which may be iden- 
tified withe the self-diffusion coefficient of a polymer. As 
dictated by the discussion in section IV later in the pa- 
per, where by comparison with data on various polymers, 
we find that a » a + a\ M n + a 2 M 2 , M n being the 
average molecular weight of the system. The constants 
in this expression are such that a(M n ) decreases as M n 
increases, as discussed above. Let us now make some 
more definitions, viz., N is the average chain length, 
the average chain density is c = pN avoga d ro /(poN), N e 
is the average chain length between consecutive points 
of entanglement, the entangled chain number density 
c e = pN avoga( i ro / (p N e ), and the monomer number den- 
sity c = pN avogadro /p Q , with p being the mass density 
of polymer, N avoga d r o the Avogadro's number, and po 
the molecular weight of the monomer. I is the length 
scale in our theory and we shall take it to be I — Ac -1 / 3 , 
where A is a parameter taken to be (2/3) 1 / 3 as it leads 
to an expression for the tensile modulus in the Gaussian 
approximation which agrees with the standard Wall the- 
ory result in the limit of short (entangled) states. We 
shall use the length I to scale all other lengths in the 
system. These ydfias are slightly similar to Stillinger's in 
another contextcjEj. By choosing our length scale in this 
fashion, it allows us to see how higher order corrections 
beyond the Gaussian approximation lead to an enhanced 
elastic modulus. In this manner we have attempted to 
relate our theory in an intimate fashion to the notion of 
entanglement. 



III. BEYOND THE GAUSSIAN 
APPROXIMATION. 



We shall now use diagrammatic methods to go beyond 
the Gaussian approximation to the structure function de- 
scribed at the end of the previous section. The reason is 
to be able to describe the crossover to a state of entan- 
glement. As discussed in the previous section, the onset 
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of entanglement is not a phase transition, but simply a 
crossover. 

Figures la and lb show the basic vertices in our theory. 
The figure captions describe the Feynman rules which go 
with these vertices. We shall compute only the first non- 
vanishing terms which arise from each of these vertices. 
The first order contribution of the cubic term is zero, as 
follows from symmetry considerations. We have to go to 
the second order in the cubic term to obtain a tadpole 
diagram which is non- vanishing, as shown in Figure 2a. 
It serves to renormalize the correlation function in the 
long wavelength limit. Figure 2b is the other setting sun 
diagram which comes from the second order contribution 
of the cubic term. It may be expanded in powers of its 
argument k. The term proportional to k 2 helps to renor- 
malize the diffusion constant g, and serves to diminish 
it, as one would expect entanglement to. Figure 3 shows 
the conventional bubble diagram coming from first order 
perturbation theory with the quartic term. It serves to 
renormalize the correlation function in the long wave- 
length limit. 

In order to render the integrals in our theory finite in 
three dimensions, we shall use the following regulariza- 
tion scheme. We shall perform an expansion of the de- 
nominator of the Gaussian structure factor in powers of k. 
We shall retain terms up to 0(k 6 ). This is essentially an 
expansion in inverse powers of a. This expansion yields 
the requisite higher order terms in the denominators of 
the Green's function to render our integrals finite, while 
ensuring that So(k) > 0. This method has the advantage 
of retaining the correct long-wavelength behavior, at the 
expense of high momentum behavior. This is acceptable, 
since we do not expect our theory to be correct at small 
wavelengths in any event. We shall shortly compare our 
method with the conventional method of renormalization 
via counter-terms. Our single-particle Green's function 
in the Gaussian approximation is now taken to be: 

S ° {k) = 1 + 2a'k 2 - 2k* + (10) 

With this definition, the net contribution from dia- 
grams shown in Figures 2a and 3 is: 



S 2a +3(a) 



I ) 5 (0) 



(11) 



Figure 2b yields a fc-dependent contribution to the self 
energy: 



s 26 (fc)= 



d 3 k' 



where: 



Sa(a) 
Sg(a) 



47 J (2tt) 3 
da + Sgk 2 + 0(k 4 ) 

1 

128 2 1 /4 7ra 3/2 
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S (k')S (\k' - k\) 



(12) 



V 2567^/2 J V2 3 / 4 a 2 



V2 



(13) 



where the integrals were performed by approximating 
the denominator of the Gaussian Green's function by 
terms up to 0(k 2 ), as this suffices to guarantee conver- 
gence of the integrals, so that there is no sensitivity to 
the higher order terms neglected. We find that the con- 
tribution from Eqn.[ll] is extremely small compared to Sa 
from Eqnjl^. This is basically what happens in the usual 
renormalization scheme, where one eliminates terms such 
as So(0) (when Eqn.^Jis used to perform the calculations) 
using appropriate counter-terms in the energy functional. 
In fact, in this scheme, So(k) decays quadratic ally with 
k, and the considerations used to obtain Eqn.n3 auto- 
matically obtain. 



With these expressions, we see that the renormalized 
value gn = g — Sg of the diffusion constant decreases 
as a is decreased. As discussed in section II, the non- 
local attractive term in our free energy was identified 
with the formation of knotty configurations, or clusterslii 
We will therefore identify gn with the dynamics of such 
clusters. Cluster dynamics have been recently observed 
using state-of-the-art techniques by Stepanek et a£j. gn 
is zero near a = 0.18. Note that a decreases as we in- 
crease N the average chain length. We thus see that as 
entanglement increases, the effective diffusion constant 
decreases, analogous to critical slowing down. We are 
not aware of explicit experimental observations regarding 
the dynamics of clusters, with which gn has been identi- 
fied, near the cross-over threshold. This effect k-pmilar 
to the considerations of Kavallis and NoolanditHl, Kroy 
and FreyEj, and Broderix et alJ^J, who study the vulcan- 
ization transition in the mean field approximation, and 
find the diffusion constant going to zero as the vulcaniza- 
tion transition is approached. They point out that this 
result agrees with experimental results. Given the anal- 
ogy between vulcanization (chemical cross-linking) and 
physical entanglement, we believe this result should be 
experimentally observable in entangled systems as well. 
Broderix et al obtained g R — > linearly with the average 
concentration. We have obtained a more complicated de- 
pendence on the concentration. We find that a as 0.18 
when the renormalized diffusion constant goes to zero. 
One could estimate a using results from the next section 
on the tensile modulus of polymers and experimental val- 
ues for polymeric elastic moduli, and then obtaining a 
value for the critical N* at which the effective diffusion 
constant goes to zero. 



The origin of 5g > can be traced back to the nonlocal 
attractive term in U e s, defined in Eqn^. This nonlocal 
attractive term, which we interpreted as giving rise to 
entanglement, is responsible for a physical signature of 
the onset of entanglement, with g^ — ► 0. 
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IV. TENSILE MODULUS. 



Y G = ZkT~c 



(17) 



It is well-known that the Wall theory result for the ten- 
sile modulus, while yielding the correct trend, does not 
agree with experimental data on moduli by a large fac-i 
tor. Edwards' application of deGenne's reptation modelEl 
provides an enhancement factor over Kuhn's result, and 
shows conceptually how entanglement leads to an in- 
crease in the stiffness of the homopolymer system. We 
will show in this section how to obtain a similar result 
in our continuum treatment. More importantly, we will 
show we can go further, and describe a crossover, as the 
mean chain length between entanglements is decreased, 
to a regime where the tensile modulus, instead of remain- 
ing fairly constant, begins to increase extremely rapidly 
as a function of decreasing N e . The reptation model is 
unable to accomplish this,Q as it assumes that the system 
is already in the entangled state, and does not account for 
inter-chain interactions, beyond assuming a preformed 
tube. 

The HelmholtzJxee energy in the Gaussian approxima- 
tion is given by:E3cj 

T G = -^kTVcS a {k = 0) = -^kTVc J d 3 xS (x) (14) 

We may represent a strained state of the system by 
the transformation x — > x! = x + u(x) in the above 
equation. This is possible because S(x) in the above 
equation represents the density-density correlation func- 
tion < c{r)c{r — x) >, so that when the system is 
strained, c(r — x) in our theory shifts to c{r l — x'), where 
f 1 = ? + u(f) . For the case of homogeneous deforma- 
tion, we shall take u(x) = e ■ a? , where e is the strain 
tensor. The strain is assumed to be volume-preserving, 
so that d 3 x = d 3 x' . Our approach is similar to that of 
Castillo and GoldbartM It is now easy to show, using a 
Taylor series expansion that the change in the pressure 



P 



(dF/dV) 



T,N 



AF G « -kTcS (k = 0)e al3 e 7 s'D a p 1 s + 0(e' 



(15) 



where the subscript G denotes the Gaussian approxi- 
mation. As needed, we can consider the expansion of Fpe 
energy to include higher orders of the strain tensorc3. 
APq is a measure of the change per unit volume of the 
energy of the system under strain. In analogy with a 
simple harmonic oscillator, the force tensor which con- 
strains the system from undergoing a strain e* is given 

by —3kTcSo(k = Q)*e : V. Thus, the stress c? required to 
produce this strain is: 

(7 = 3kTcSo{k — 0)*e : 2? (16) 

One may now readily write down the tensile modulus 
as: 



This is identical to the well-known Wall theory result 
obtained by in the limit of short (uncntangled) chains. 
And it's origin is purely entropic. Our goal is to go be- 
yond this result, and to do that, we shall dress the bare 
propagator So(k) using the diagrams shown in Figures 2 
and 3. This immediately leads to the renormalized result 
Y R : 



Y R (a)= 3kTcS R (k = 0) 



3kTd 



1 - £ 2a +3(a0 - Sa{a) 



(18) 



The first of these equations is similar to the connection 
made between the structure factor in the lpog wavelength 
limit and the bulk modulus by KirkwoodEj. 

The result of plotting the entanglement factor Z = 
Yr/Yq — 1 as a function of a(N) is presented in Figure 
4. We see that the enhancement factor, which is fairly 
constant above a — 0.2008, begins to increase dramat- 
ically below this value of a = 0.2008. As discussed in 
section II, decreasing a is equivalent to increasing TV, the 
average number of links in a polymer. And increasing N 
is associated with increasing entanglement. 

To see better the connection between our approach and 
the underlying chain parameters, we identify the pref ac- 
tor (Z(a) + 1) multiplying the Wall theory result with 
the enhancement obtained within the reptation modelu, 
i.e., Z(a) + 1 = (A r fe 2 /a 2 ), where b is the monomer 
length, and a is the diameter of the tube in the rep- 
tation model. This identification gives us a relation be- 
tween a and the parameters of reptation theory. It also 
gives us a relation between our theory and the under- 
lying chain parameters such as N, the average number 
of links, and N e , the average number of links between 
successive points of entanglement, viz., A*"" 1 = (b 2 /a 2 ). 
In fact, we can provide such an analytic, approximate 
relation in the following manner. Based on numerical es- 
timates of 6>o(0), we find that in the the highly entangled 
state, |(5a(a)| >> |S2 a +3(a)|- Then, using Eqnjl^, we 
immediately obtain: 



a(N e , N)fa k(1 — N e /N)~ 2 ^ 3 



1 



128 2V4 n 



2/3 



(19) 



It should be noted that this particular scaling rela- 
tion Eqn.[l9] holds only in the limit of highly entangled 
states. While we have denoted the explicit dependence 
on N e and N separately, it should be noted that N e it- 
self depends on the chain length N. More generally, the 
advantage of our theory is that the enhancement factor 
[Z(a(N)) + 1] can change continuously from a value of 
unity in the unentangled state, to a fairly large number 
as the system becomes increasingly entangled. As dis- 
cussed earlier, the reptation model on the other hand 
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presupposes the formation of tube-constraint caused by 
polymer chains surrounding any given chain. As such it 
applies only in the highly entangled state. 

We introduced our theory in section II in a phcnomeno- 
logical fashion, and so the best way to establish the va- 
lidity of our theory is to compare the results of our the- 
ory with experiment, the connection with the underlying 
chain parameters (Eqn.|l9|) notwithstanding. However, 
we are not aware of experimental data on the elastic 
properties of entangled polymers detailing in a precise 
manner the cross-over to the entangled state. 

Nevertheless, we-were able to find an experimental pa- 
per by Onogi et alEZl which provides systematic viscoelas- 
tic data on amorphous polystyrene in the melt state, over 
a wide range of frequencies and molecular weights (and 
having a low degree of polydispersity). We were able to 
obtain static shear moduli from the plateau moduli given 
in this paper. The conventional assumption of incom- 
pressibility then says that Young's modulus is three times 
the shear modulus. Armed with these data, we were 
able to obtain a parameterization of a as function of the 
molecular weight M n (= noN) (see Table 1 and Figure 5). 
This parameterization is consistent with the stress-strain 
data provided by Bicerano et aH. This parameterization 
is also consistent with Eqn.[l9| Equation [l9] implies that 
for very large molecular weights, a — > k ~ 0.01635. From 
Figure 5, we see such values for a occurring at the high 
end of the range of molecular weights. The parameteri- 
zation depicted in Figure 5 may be thought of as a power 
series expansion of Eqn. 

The cross-over between the un-entangled state and 
entangled state appears to take place between M n w 
113,000 and M n w 60,000, as the plateau in the stor- 
age moduli measurements disappears somewhere between 
these two values of the molecular weightsEH. Unfortu- 
nately, this cross-over is not described well by these data, 
as measurements in the cross-over region are unavailable. 
This could well be due to the critical slowing down which 
our theory predicts will occur during the cross-over (see 
section III). As can be seen from Table 1 and Figure 5, 
the data develops noise at the low end of the range of 
molecular weights. And this is consistent with the pre- 
diction of critical slowing down, as the system takes very 
long to settle into a metastable state. Since we have a 
cross-over phenomenon, it is difficult to pinpoint a sin- 
gle value of M n at which the transition to the entangled 
state occurs. But based on Figure 6, in which the log- 
arithm of the enhancement factor has been plotted as a 
function of M n , we see that it starts to become appre- 
ciable (>> 1) in the neighborhood of M n ~ 50, 000. The 
theory thus predicts that it is in this neighborhood that 
the cross-over takes place. It would be very interesting 
to see more experiments performed in the future, say for 
amorphous polystyrene, with M n in the range just dis- 
cussed, to see if the path of the cross-over agrees with 
that shown in Figure 6. It is even more important to do 
such experiments to ascertain if the critical slowing down 
predicted in section III does indeed occur. 



V. CONCLUSIONS 

We postulated an extension of the Cahn-Hilliard func- 
tional to describe entanglement in polymers. We ex- 
tended the Cahn-Hilliard functional with two terms. One 
is an attractive nonlocal term which describes the effect 
of entanglement, and the other a local repulsive term in- 
dicative of excluded volume interactions. We showed in 
Appendix A how this extended functional can be derived 
using notions from gauge theory. Using field theoretic 
techniques to go beyond the Gaussian approximation, we 
showed that the onset of entanglement is a crossover phe- 
nomenon, signaled by the effective diffusion constant go- 
ing to zero. We developed a simple model to connect the 
single parameter in our theory with the parameters of the 
underlying chains by comparison with available experi- 
mental data on amorphous polystyrene melts. A reason- 
able estimate for the critical partial chain concentration 
at which this crossover occurs was obtained. While this 
is consistent with available data, further experiments to 
determine the details of this cross-over were suggested, 
especially to ascertain if the critical slowing down pre- 
dicted in this paper does indeed occur. 
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APPENDIX A. 

In this appendix, we shall show how to derive Eqn^, 
along with the choice of parameters made in section II, 
using notions from gauge theory. Let us start with: 



Puo(c(b)) = (I) 



dc(s) dc(s) 
dsi ds. 



(A.l) 



where the variables have been defined in section II. We 
will use this form as our starting point to generate a more 
complete energy functional using gauge invariance. 

From Eqn. A.l we see that uq is invariant under global 
translations of c(s), i.e. under c(s) — > c(s) + h where h is 
a constant. And the appropriate group to consider is Ti. 
The physical origin of this group can be traced back to 
the fact that the quadratic (positive, semi-definite) form 
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of the energy density is dictated by expanding the inter- 
nal energy around a minimum, in a Landau-like fashion. 
Physicality of Ti transformations demands that c' + c e > 
0, where d denotes the deviation of the specie concentra- 
tion from its average, and that number conservation is 
guaranteed. These physical constraints will be incorpo- 
rated into the evaluation of the partition function, after 
gauging Uq, in a manner similar to applying gauge con- 
straints in QFT. 

Our physical motivation for seeking local gauge invari- 
ance under Tj is the same as that of Yang and MillsE], 
and in quantum electrodynamics (QED), where one ob- 
serves the invariance of the noninteracting Lagrangian 
under certain global transformations. One then demands 
covariance of the theory when these symmetry operations 
are local i.e., when the transformations are space-time de- 
pendent. A reason for this, as given by Yang and Mills, 
is that one can now freely interchange between the fields 
as one moves through space and time, while leaving the 
physics covariant. It is important to note that gauge the- 
ory in QFT is not a result of the fact that the phase of 
the field is not measurable. In fact, Aharonov and Bohm 
showed many years ago that the phase in quantum me- 
chanics is indeed observable. Local transformations un- 
der Ti generate concentration fluctuations which arise 
from entanglement. ■— ■ 

Following Yang and MiUs,E3 local gauge invariance of 
wo under Ti motivates us to define new fields b, which 
have invariance properties appropriate to Ti. We define 
a covariant derivative — > (^- + gr&j), where r = 
is the generator of Ti, q is a 'charge', or equivalently, 
a coupling constant, and the 6-fields are analogs of the 
magnetic vector potential in electrodynamics. In our pre- 
vious theory of self-assemblyEJ, gauge fields arose from 
the underlying covalent bonds between the two species 
in the system. In the present case, where we wish to de- 
scribe entanglement, the gauge fields are to be thought 
of as arising purely from statistical considerations alone. 
On the other hand, chemical cross-linking would provide 
a physical origin for the Wields in vulcanized systems. 
The energy functional for the 6-fields is defined a la Yang 
and Mills, via the minimal prescription. With this, our 
original internal energy density is transformed into: 



J(c) = gqVc 



/3u — ► (3u = (3u + (3u int + Puym 



(A.2) 



where Ui n t refers to the interaction energy density, and 
uym is the energy density associated with the Yang-Mills 
fe-ficlds alone. Equivalently, we may define the total en- 
ergy functionals associated with these energy densities: 



f3U ^(3U = PU + f3Uint + PU YM , 



where 



font = J(c) ■ 6(S) + / b( S ) • 6(8) 



(A.3) 



/ = 



2 



We need one more definition for completeness: 



Puym 



dSj 



dsj 



dsi 



(A.4) 



(A.5) 
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(A.6) 



This equation can be cast into the following form: 



; "> w = - ( g 



(A.7) 



EqnA.7 is obtained via an integration by parts, in the 
transverse gauge. Since we are dealing with an Abelian 
gauge theory, it is permissible to insert this transverse 
gauge manually, without resorting to the formal machin- 
ery of Faddeev and Popov. 

Again using the transverse gauge and integrating by 
parts, it is clear that: 



d 3 sVc(s) • b{s)= - / d 3 s c(a) [V • b(s)] = 



i / d 3 sVc(s) • b(s) 



with 



(A.8) 



It is this crucial identity which allows us to get the 
precise form for Eqn.^j, which we motivated in section 
II using an intuitive approach. It is the nature of the 
Ti group which permits this manipulation to go through 
successfully. In our earlier papertL3, where we used the 
SO (2) group, such a manipulation would not have availed 
us any advantage. 

Note that we are utilizing a non-relativistic version of 
the Yang-Mills procedure, since we are only concerned 
with time-independent problems. Furthermore, since we 
are concerned with translations in Ti, there is only a sin- 
gle generator to contend with, so that the resulting func- 
tional is only quadratic and not quartic in the 6-fields. 

It is important to emphasize that the usual application 
of the Yang-Mills procedure in QFT implies the existence 
of fundamental interactions. In our case, we are applying 
the principle of local gauge invariance at the mesoscale. 
Consequently, we do not expect to discover any new fun- 
damental interactions by using gauge invariance. Rather, 
we interpret the new Wields as yielding correlations be- 
tween the concentration fields. These correlations could 
also be thought of as effective interactions, which arise 
at the mesoscale from the underlying electrostatic inter- 
actions between molecules. 

The partition function we need to evaluate is now: 

Q'= fvc9 (c) ]J Vb k exp - 0(U O + U mt + U YM ) 

(A.9) 
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Equation |A.9| is a functional integral, where the step 
functions denoted by 9 imply that we must restrict inte- 
gration to positive semi-definite values of the fields. 

Since the 6-fields appear only quadratically in the 
above functional, it is straightforward to integrate over 
them, and obtain an effective int erna l energy functional 
involving only c, upon using Eqn. A. 8. The result is: 



= pU + 



1 
4 



■0AU ( 



err 



d J s / dV Ji(c(s)) 



J,(c(s')) 



Note that in doing so, we have ignored an overall triv- 
ial normalization constant that appears in the evaluation 
of the partition function Q'. This is permissible, as this 
factor cancels during the evaluation of averages of ob- 
servable quantities. 

To reveal the physics in this effective functional, we 
perform some straightforward algebra to write our result 
as: 



(3U cS = (3U + 



J d 3 s c(s)c(s) 



2tt 



d J s / dV c(s) 



exp(- v /2a 2 /g\s -s'|) 



s - s' 



c(s') (A 



where a 2 = g 2 q 2 /2 (we shall use units in which q=l). 
Note that U e s is quadratic, the generator of Ti making 
sure that higher order terms do not appear in our func- 
tional. Equation A. 11 is one of the main results of our pa- 
per, and provides a deeper motivation for the model de- 
veloped in section II on intuitive grounds. As discussed in 
section II, the form of E quatio n [A. 11 guarantees number 
conservation. Equation A. 11 shows that entanglement 
may be understood in the context of a mesoscopic gauge 
theory. Note that the two terms we just discussed have 
signs opposite those of corresponding terms in theories 
of self-assemblylii We thus see that using Ti instead of 
SO (2) in the previous theoryEJ has led to a qualitatively 
different theory. Finally, this approach provides an alter- 
native gauge theory oLpolymer entanglement than the 
one given by BreretonLJ. In this paper, Brereton works 
explicitly with entangled strands of polymers and pro- 
vides a static theory which yields e.g. a renormalized 
expression for the radius of gyration. It would be inter- 
esting to see how Brereton's theory could be generalized 
to include dynamics as well. 



3 M. Doi, Introduction to Polymer Physics, Oxford Univer- 
sity Press, Oxford (1996). 

4 M. Doi, S.F. Edwards, The Theory of Polymer Dynamics, 
Oxford University Press, Oxford (1986). 

5 S. Prager, H.L. Frisch, J. Chem. Phys., 46, 1475 (1966). 

6 K. Koniaris, M. Muthukumar, J. Chem. Phys., 103, 7136 
(1995). 

7 Y. Termonia, P. Smith, Macromolecules, 20, 835 (1987); 
Y. Termonia, Macromolecules, 29, 4891 (1996). 

8 J. Bicerano, N.K. Grant, J.T. Seitz, K. Pant, to appear in 
(Aj'Steromolecules (1997). 

9 T. Holtzl, H.L. Trautenberg, D. Goritz, Phys. Rev. Lett., 
79, 2299 (1997). 

10 S.M. Chitanvis, Phys. Rev. E 
57, 1921 (1998); also available in post script format at 
http://www. tl2. lanl.gov/~shirish/preprints. html. 

11 J.W. Cahn and J.E. Halliard, J. Chem. Phys.,28,258, 
(1958). 

12 T.A. Kavallis and J. Noolandi, Phys. Rev. Lett., 59, 2674 
(1987). 

13 K. Kroy and E. Frey, Phys. Rev. Lett., 77, 306 (1996). 

14 H.E. Castillo and P.M. Goldbart, preprint, cond- 
mat/9712050. 

15 R.T. Deam and S.F. Edwards, Phil. Trans. R. Soc, 280A, 
317 (1976). 

lf J P.G. deGennes, Scaling concepts in polymer physics, (Cor- 
H) nell University Press, Ithaca, N.Y. 1979). 

17 H.B. Callen, Thermodynamics, p. 106, J. Wiley & Sons, 
Inc., N.Y. (1960). 

18 H.J. Raveche, J. Chem. Phys. 55 2242, (1971). 

19 P. Stepanek and Wyn Brown, Macromolecules 31, 1889 
(1998). 

20 F.H. Stillinger, J. Chem. Phys., 78, 4654 (1983). 

21 C. Carraro, Physica A, 236, 130 (1997). 

22 K. Borderix, P.M. Goldbart, A. Zippelius, Phys. Rev. Lett., 
79, 3688 (1997). 

23 P. Ramond, Field Theory: A Modern Primer, page 126, 
The Benjamin/Cummings Pub. Co., (1981). 

24 J.J. Binney, N.J. Dowrick, A.J. Fisher, M.E.J. Newman, 
The theory of critical phenomena: An introduction to the 
renormalization group, Oxford Science publications, Ox- 
ford (1995). 

25 N. Triantafyllidis and S. Bardenhagen, Journal of Mech. 
and Phys. of solids, 44, 1891 (1996). 

26 J.G. Kirkwood, R.J. Goldberg, J. Chem. Phys., 18, 54 
(1950). 

27 S. Onogi, T. Masuda and K. Kitagawa, Macromolecules 3, 
109 (1970). 

28 W.M. Graessley, The entanglement concept in Polymer 
Rheology, Adv. Poly. Sci. 16, Springer- Verlag, Berlin 
(1974). 

29 C.N. Yang and R.L. Mills, Phys. Rev., 96, 191 (1954). 

30 M.G. Brereton, J. Molec. Struct. 336, 191 (1995). 



1 C.G. Moore and WF. Watson, J. Polymer Sci., 19, 237 
(1956). 

2 P.J. Flory, N. Rabjohn, M.C. Shaffer, P. Polymer Sci., 4, 
225 (1949). 



8 



TABLE I. Fitted values of a obtained in comparison with experimental data from Onogi et al 
on amorphous polystyrene melts at 160 °C. Young's modulus is denoted by Y. 



Molecular weight 



Yexptl (psi) 



a(fitted) 



581,000 
513,000 
351,000 
275,000 
215,000 
167,000 
113,000 



93 
84 
66 
65 
62 
56 
106 



0.0166 
0.0167 
0.0171 
0.0173 
0.0177 
0.0184 
0.0179 



FIG. 1. (a) is a pictorial representation of the cubic term in A. Each leg corresponds to a factor 
of c, the field. The intersection of the three legs symbolizes a factor of 7 = 1/6, the coupling 
constant, (b) is a pictorial representation of the quartic term in A. A factor of —1/12 is to be 
inserted at the intersection. 

FIG. 2. (a) represents the tadpole diagram which is crucial in our calculations, (b) represents 
the setting sun diagram. Both (a) and (b) are second order contributions to the correlation function 
coming from the cubic interaction term, the first order corrections being null. 

FIG. 3. This figure represents 1-loop (bubble) contribution from the quartic interaction term 
in A. 



FIG. 4. This is a plot of Z = Yr/Yg — 1 as a function of a. Notice that the factor is virtually 
constant above a — 0.20, followed by a dramatic increase below this value of a. Approximations 
employed in the calculation cause the entanglement factor to diverge at a m 0.01. Remember, 
decreasing a corresponds to increasing entanglement. 

FIG. 5. This dots are the values of a required to achieve agreement with experimental values of 
moduli of amorphous polystyrene melts for a range of M n . The solid line indicates a least squares fit 
obtained using a(M„) = a +a l M n +a 2 Mi, with a = 0.0189, 01 = -6.18xl0~ 9 , a 2 = 3.89xl0~ 15 . 
Caution should be exercised in using this expression to estimate a too far from the range in which 
the fit was made. 

FIG. 6. Plot of Log[.2] as a function of the molecular weight M n . Notice that Z begins to attain 
nontrivial values above M n ~ 5 x 10 4 . The curve indicates the path of the cross-over from an 
un-entangled state to an entangled state. 
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Chitanvis, Fig. [l] 





(a) 



(b) 



CM 



12 



Chitanvis, Fig. [| 
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0.0185 



0.0180 



0.0175 



0.0170 




n rn 



